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We study the ofT-equihbrium effects of inflaton on the dynamics of primordial perturbations in the 
0{N) model. A self-consistent off-equilibrium formalism is employed to investigate the evolution 
of the inflationary background field and its fluctuations with the back-reaction effects. We find 
two observable remains left behind the off-equilibrium processes: the running spectral index of 
primordial density perturbations and the correlations between perturbation modes in phase space, 
which would serve as the imprints to probe the epoch of inflation, even beyond. 
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In the inflationary scenario, the primordial pertur- 
bations of the universe originate from vacuum fluctua- 
tions of the scalar fleld(s), the inflaton 4>, driving the 
inflation. If the dynamics of the fluctuations is ap- 
proximated by single massless free field during the in- 
flation, the power spectrum of curvature perturbations 
for the mode k, VTi{k) is roughly fe-independent 
However the fc-dependence, or the running of the spec- 
tral index, revealed by the recently released Wilkin- 
son Microwave Anisotropy Probe (WMAP) data can 
be as large as dn{k)/d\iik = {d/d\nk)'^ InV-jiik) — 
-0.055 to - 0.077 la. Hence, the origin of the pri- 
mordial perturbations cannot be solely accommodated 
with the quantum fluctuations of single free field. In 
a slow-roll inflation, one has d/d\nk ~ {l/H)d/dt = 
-{l/87rG){V'{(t>)/V {(/))) {d/d(j)), and thus the derivative 
of the power spectrum (d/dlnk)"' InVnik) will no longer 
be negligible for the second order (n = 2) as long as 
the third or higher derivatives of the inflaton potential 
becomes substantial. To fit in with the running of the 
spectral index, models beyond single scalar field with 
quadratic potential ¥{(/)) have been proposed accord- 
ingly 1^. In this context, the running index of the per- 
turbation power spectrum is considered as an essential 
feature due to the interactions or the self-interactions of 
the infiaton(s). 

It should be pointed out here that in deriving the 
derivative of the power spectrum above, the effective po- 
tential 1^(0) is involved. However, this approach may be 
problematic, when in particular, the one-loop effective 
potential becomes complex where the background field (j) 
is in the region of V {(j)) < |3|. The imaginary part 
of the effective potential would inevitably lead to a dy- 
namically unstable state . This so-called "spinodal in- 
stability" will allow long wavelength fluctuations to grow 
non-perturbatively. Therefore, the primordial perturba- 
tions must be dealt consistently with a different method 
to account for the ampliflcation of vacuum fluctuations in 



the presence of spinodal instabilities. This motivates us 
to study the self-interaction effects of the inflaton on the 
dynamics of primordial perturbations within a context of 
the self-consistent off-equilibrium formalism. We will use 
the 0{N) model of inflation as an example. The 0{N) 
model with spontaneously broken symmetry has been ex- 
tensively used in modelling the quantum off-equilibrium 
processes in the early universe, as well as the chiral phase 
transition in relativistic heavy ion collision 0, H, ■ 
Here, we will focus on searching for the possible observ- 
able imprints caused by the off-equilibrium evolution of 
inflationary primordial perturbations jlfl] . 

Consider the dynamics of an inflation driven by a field 
$ of the 0{N) vector model with spontaneous symmetry 
breaking. The action is defined by 
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where • ^) is a self-interaction potential given by 
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As an inflaton, the components of the field generally 
are represented as $ = (cr, 7?) , where 7? represent A^ — 1 
scalar fields. The cosmic inflation is characterized by 
the state in which the component a has a spatially ho- 
mogeneous expectation, i.e. a can be decomposed into 
a background held plus fluctuations around the back- 
ground held as 



(3) 



with the expectation value ((t(x, t)) = y N(j){t) and thus, 
(x(x, t)) = 0. During the inflationary epoch, the back- 
ground space-time can be described by a spatially fiat 
Fricdmann- Robertson- Walker metric, 



cr(x,i) = ^N(j){t) +x{^,t) 



ds^ = g^^dx'^dx'^ — dt^ — a^[t)5ijdx^dx-' , 



(4) 
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where the scale factor a = exp{Ht) with the expansion 
rate H = ^SnGp/S determined by the mean of energy 
density of the inflaton field, p. 

In order to take account of the growth of fluctuations 
due to spinodal instabilities as we will see later, we will 
employ the method of the Hartree factorization, which 
approximates the potential V with an effective quadratic 
potential while keeping N finite 7,9]. Then the equation 
of motion of the background field (p can be directly ob- 
tained from the Hartrce-factorized Lagrangian by means 
of the tadpole condition (x(x, t)) = given by 

4>it) + 3H<p{t) + [Ml (t) - (t)] = . (5) 

We then decompose xi^) ^^id 7f(t) in the Fourier basis. 
In the Heisenberg picture, one has 



0.6 - v{A/2) 



(6) 



where Oik, &k and a|j^, are the creation and annihi- 
lation operators which obey the commutation relations. 
The equations of the mode functions fx.kit) ^^^d /ir,k(0 
can be found from the Heisenberg field equations as fol- 
lows: 



—2+^H- + -+Ml{t) 



dt 



dt 



d ^ d k , , , 



fxMit) = 0, 



UM{t) = 0. (7) 



The time dependent mass terms are give by 



(8) 



where {ip'^) is defined by (t?^) = {N — The ini- 

tial conditions of the mode functions can be specified as 
effective free massive scalar fields in an expanding uni- 
verse 0: 



/x,k(0) = 

A,k(o) - 



1 



2[F + M2(0)] 
1 



■ 4k(o) = -* 



/.,k(0) = -i 



+ Af2(0) 



fc2 + M2(0) 



where we have set a(0) = 0. The values of M2(0) and 
(0) depends on the details of the onset of the inflation. 
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FIG. 1: The evolution of J^^, 

V 2m' 

spectively vs. t {in units of m^^ ) for initial conditions (j}{0) 
0, 4>{0) = 0,and the parameters h = S. ^ 2.0, A = 10~" in 
the case of A'' = 4. 



Finally, to close these equations self-consistently, the 
terms (x^)(0 {ilP'){t) in the mass-squared [Eq.®] 
can be determined by the mode functions which involve 
the divergences associated with the loop integrals. To 
get the self-consistent renormalized quantities defined as 

{x')R(t) = rsi/x.kWp ' 
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one has to absorb the divergences into renormalization 
of the bare mass and the bare coupling constant. How- 
ever, for a weak coupling A < 10"^'* in a typical inflation 
model, the logarithmic divergences can be neglected. 

Figure 1 plots a self-consistent solution of the inflation- 
ary background field (j){t) and the fluctuations {tp'^){t), 
{x^){t), in which the parameters are taken io he h — 
H/m = 2.0, A = 10-1'*, TV = 4, and M2(0) = M2(0) ~ 
rri^ > 0. In fact, the value we choose for H/m above 
can be determined by the energy density of the infla- 
ton field. During the inflation, the energy density of 
the inflaton is dominated by that of the potential en- 
ergy that leads to {H/mf « Nrn'^/M^^X w 0(1) for 
N « A « 10"", and the inflaton mass m « 10*2 

Gev. We see that ~ and (V'^)W ^ (x^)(i) - 
during the first stage of the inflation as i < AOm~^ in 
the spinodal regime where the initially positive M^. and 
M2 turn negative. The spinodal instabilities become im- 
portant eventually, and leads to the significant growth of 
both {ip^)it) and (x^)(t) starting at t = 50m-^ Then, 
the spinodal condition is only weakly satisfied over a 
span of 70m~^ < t < 80m~^ where the background cj) 
starts the rapid falling into the valley of the inflaton po- 
tential. Finally, the process of spinodal instability ter- 
minates at time given by M^{te) ~ 0. It renders 
te — (70 — 80)m~i. The background inflaton cf) sets in 
the minimum of a broken symmetry phase at the time to 
where (j) — •\/2/Am, which signifies the end of inflation. 
Therefore, the number of total inflationary e-foldings is 
iVe w Ht = 2mt ~ 160. It is evident that most of energy 
density of the background field is transferred to the 
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fluctuations {ip"^) via the production of tt modes during 
the inflation. Since M^{te) — 0, these tt modes are in 
fact the massless Goldstone modes of the broken symme- 
try phase. The Goldstone theorem is fulfilled dynami- 
cally. Although the Hubble parameter H is treated here 
as a constant, the energy density of the inflaton evolves 
dynamically. In fact, the time dependent H can be de- 
termined by the evolution of the inflaton fleld as the 
effects of quantum fluctuations (« 10~^ ) on the dynam- 
ics of the Hubble parameter is negligible. The approxi- 
mation of the constant H can be argued to beak down 
when the (p{t) reaches the maximum value at < = 70m^^. 
A truly self-consistent approach has to involve the cor- 
rect dynamics of the Hubble parameter through the semi- 
classical Einstein equations as in Ref. "t^. However, the 
qualitative features of the above solutions are shown to 
remain the same. 

Using the numerical solutions, we compute the power 
spectrum of primordial perturbations. The power spec- 
trum of primordial perturbations is described by V{k) — 
(l^kP) where the mass density perturbations are deter- 
mined by the gauge invariant quantity ^3 



4 = 



Sp 



P + P 



(10) 



k=aH 



The mass density fluctuation Sp originates from the fleld 
fluctuations tt^ and x given by 
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The term summing up the energy density and the pres- 
sure, p + p, in Eq. ifTUI) can be obtained by 
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(12) 



The spectral index n{k) — 1 = dP{k)/d\nk and its k- 
dependence dn{k)/d\ak are computed numerically as 
shown in Fig. 2. The value of n{k) varies from unity at the 
larger scales as i < 60m~^ to about n = l-t-4(j/ — 3/2) ~ 
1.32 at the smaller scale as i > 80to~^. Apparently, there 
is a significant index-running dn(k) / d\nk ~ 0.015 in the 
wavelength range corresponding to the horizon-crossing 
times during 60 < mt < 80. The running of the spec- 
tral index is in fact due to the energy transfer from the 
inflationary background fleld to the fluctuations which 
certainly can not be obtained from the classical effective 
potential approach. 




FIG. 2; Spectral index n{k) and its running dn(k) / d\nk vs. 
Infc ( fc in units of m ) with initial conditions (j!)(0) ~ 0, <^(0) = 
0, and the parameters h = — = 2.0, A = 10~^* in the case of 
TV = 4. 



The correlation of fluctuations caused by the off- 
equilibrium process is generic, because the evolution of 
fluctuations is typically with respect to a time-dependent 
background. When perturbations of scale k cross the 
horizon at a time i, the equal-time two point correlation 
function between (x, t)-(x',t) will yield the correlation 
between two space-scale modes (x, fc)-(x', k) via the map- 
ping formula k = a{t)H in Ea. (|10|l . The discrete wavelet 
transform (DWT) is designed to do such space-scale [x-k] 
decomposition p!3 |. 

In the formulation of DWT, there are two sets of spa- 
tially localized bases given by the scaling functions |j, 1)^, 
and the wavelet functions |j,l)M>; both are characterized 
by the indices j and 1 where j = 0, 1,2, ... stands for a 
scale from kj to kj + Akj in which kj = 2-k'2? jL and 
Akj = 2-112^/1. The index 1 = 0, l,...,2j - 1 denotes 
the location of the spatial point within L\/2^ < xi < 
L(l + l)/2j. These bases are complete, and they sat- 
isfy the orthogonal relations s(j', l)s = <5jj'(5i,i', and 
^ (j', l'|j, 1)^ = '5jJ"^i,i'- Then , the two-point correlation 
functions of the density contrast 5 can be rewritten in 
terms of the DWT bases as 



d^k 



^(<5(k,i)J(-k,i))V'ja(k)^j;,,(k)^ 

d^k . „ 

{5M - / ^(<5(k,t)j(-k,t))0j,i(k)0;,,(k). 



(13) 



where V'jj(k) = (k|j,l)s,and 0jj(k) = (k|j,l)„ are the 
scaling functions and the wavelet functions in the k- 
representation. The time t is taken to be tj specified 
by the relation 2tt2^ / L = k = aH . Since a ~ exp{Ht), 
and one has 
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Thus, Ea. (|13l) can be used to determine the correlations 
between fluctuations at different spatial points 1 and 1', 
both crossing out of the Hubble horizon at the same time 
tj during the inflationary epoch. 



4 





•— • j=S. AT^ > 
*—• j=6. AT^ > 




0.2 (l.-l 0(j {)S 




•— • j=4. Af=0 
j=5. M^,=0 
•— • j=6. M^,=0 





„i , \ , \ — — 1 

I) 0.2 0.4 0.6 0.8 1 

r=L(Wl/2' 



FIG. 3: The normalized mode-mode correlation functions of 
density fields in phase space are plotted against the separation 
r ( in units of ) of two perturbation modes with respect 
to various scales j under two difi'erent sets of initial conditions. 
The parameter L is taken to be 2m~^. 



As a numerical example, the normalized two-point 
space-scale correlations of density perturbations defined 
as (<5j,ii5j,i')/('5|i) under different initial conditions are 
plotted in Fig. |31 in which the parameters are taken to 
he H ~ 2m and L — 2m^^. From Ea. ((n|l . one has 
tj ~ (1.63 — 2.33)m^^ for j = 4 — 6 which corresponds 
to the number of e-foldings Htj ~ 3.3 — 4.7. Since the 
inflation under consideration lasts from i = to about 
t — 80m^^, the correlations in Fig. 13 actually probe the 
inflaton dynamics at the beginning of the inflation. With 
M^(0) = M.^(0) ~ > 0, the space-scale correlation 
function approaches zero drastically as the distance r be- 
tween the two modes increases; while the correlation de- 
viates from zero and lasts for a long range for perturba- 
tions with = M"^ = initially. It indicates that the 
evolution of the correlation depends upon the initial con- 
ditions of the inflation. FromEq. H14l 'l. the v'-dependence 
of the mode-mode correlation can probe the evolution (t- 
dependence) of the primordial perturbations. Thus, one 



is capable of exploring the physics of the very early uni- 
verse by means of the mode-mode correlations in phase 
space. 

In summary, we study the off-equilibrium effects on 
the primordial perturbations in the 0{N) model. When 
the interaction or self-interaction of inflaton becomes im- 
portant, one must consider the evolution of the inflation- 
ary background field and its fluctuations with the back- 
reaction effects using a self-consistent off-equilibrium for- 
malism. We find two observable remains left behind the 
off-equilibrium processes which would serve as the im- 
prints to probe the epoch of inflation. The first one is the 
running spectral index of primordial perturbations. We 
find that the running spectral index depends essentially 
on the rate of the energy transfer from the background 
field to the inflaton fluctuations. The second remain is 
the correlation between phase space modes of the density 
perturbations. Under the influence of the self-interaction, 
fluctuations created from the background field arc no 
longer white noises. Moreover, since the evolution of the 
correlation depends upon the initial conditions of the in- 
flation, the mode-mode correlation of density perturba- 
tions also provides a window to study the dynamics of 
the self-interaction as well as the onset of the inflation. 
Thus, we may expect that the non-trivial mode-mode 
correlation in the phase space is detectable via a DWT 
analysis on the CMB temperature map, or other observa- 
tions on the large scale structure relevant to the density 
perturbations |l4| . 
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